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Abstract 

We present an algorithmic proof of the Cartan-Dieudonne theorem on generalized 
real scalar product spaces with arbitrary signature. We use Clifford algebras to 
compute the factorization of a given orthogonal transformation as a product of re- 
flections with respect to hyperplanes. The relationship with the Cartan-Dieudonne- 
Scherk theorem is also discussed in relation to the minimum number of reflections 
required to decompose a given orthogonal transformation. 
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1 Introduction 



The Cartan-Dieudonne theorem is a fundamental result in the theory of met- 
ric vector spaces. It states that every orthogonal transformation is the com- 
position of reflections with respect to hyperplanes. The classical proofs of the 
Cartan-Dieudonne theorem use induction on the dimension of the vector space 
and are not constructive. See [I] and [2]. Recently, Uhlig [3] presented a con- 
structive proof of the Cartan-Dieudonne theorem for the case of vector spaces 
with a positive definite inner product, and also a constructive proof of a weaker 
version of the theorem for generalized scalar product spaces of signature (p, q) 
[3j Thm.3]. 

The matrix representation of a reflection with respect to a hyperplane is called 
a Householder matrix [3]. The analogues of Householder transformations on 
spaces with a non-degenerate bilinear or sesquilinear form are studied in [3]. 
Householder matrices are also used in Gallier's book [51 Ch. 7], which discusses 
the Cartan-Dieudonne theorem for linear and afline isometries, including ap- 
plications to QR decomposition. 

In the present paper we present an alternative proof of the Cartan-Dieudonne 
theorem for generalized real scalar product spaces of arbitrary signature. The 
proof yields an algorithm for the factorization of a given orthogonal transfor- 
mation as a product of reflections with respect to hyperplanes. This work is 
a generalization to spaces of arbitrary signature of a previous one [5], where 
we provide an algorithmic proof of the Cartan-Dieudonne theorem in IR n valid 
over the fields Q, K or C. 

In the theory of Clifford algebras there is an alternative way to find the image 
of a vector under a reflection with respect to a hyperplane. This is done us- 
ing vector multiplications under the rules of the Clifford algebras. This is the 
method that we propose for the computations involved in the factorization 
algorithm. Aditionally, our approach produces an alternative way for calculat- 
ing the Householder matrices with respect to orthogonal bases, but we don't 
use these matrices in the numerical examples. 

In section 2 we present some results about vector spaces with non-degenerate 
bilinear forms. We include some results about Artinian spaces because they are 
important for the development of our proof of the Cartan-Dieudonne theorem. 
In section 3 we present the proof for the case of spaces with a symmetric non- 
degenerate bilinear form of arbitrary signature. In section 4 we propose the 
use of Clifford algebras as a computational tool to obtain the reflections that 
give the factorization of a given orthogonal transformation. In section 5 we 
present some examples of the factorization of orthogonal matrices. Finally, in 
section 6 some conclusions are given, including a comment about the relation- 
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ship with the Cartan-Dieudonne-Scherk theorem, and the minimum number 
of reflections required to decompose a given orthogonal transformation. 



2 Generalized scalar product spaces 

In this section we present some basic results concerning real vector spaces 
equipped with a non-degenerate symmetric bilinear form. We call such spaces 
(generalized) scalar product spaces. They are also known as metric vector 
spaces or real orthogonal spaces. The proofs of the results presented here can 
be found in [2]. 

Definition 2.1 Let X be a real vector space and let B : X x X — )■ K be a map 
that satisfies the conditions 

(Bl) (Bilinearity) For allv,v',w,w' G X and for all X G K 
B (Xv + v', w) = XB (v, w) + B (V, w) , 

and 

B(v, Xw + w') = XB (y, w) + B (y, w') . 

(B2) (Symmetry) For all v , w G X , B (v, w) = B (w, v). 

(B3) (N on- degeneracy) For each non-zero v in X there exists w in X such 
thatB(v,w) j£ 0. 

Then the pair (X,B) is said to be a (generalized) scalar product space. 
Definition 2.2 Let X be a real vector space and B a bilinear form on X. 

(1) The vectors u,v in X are orthogonal if B (u, v ) = 0. 

(2) A vector u in X is called isotropic if B (u,u) = 0. 

(3) We say that u is invertible if u is not isotropic, that is, if B(u,u) ^ 0. 
(This terminology will be justified in section 4-) 

(4) Let W and V be vector subspaces of X . We say that W and V are or- 
thogonal if B (v, u) = 0, for all u G V and v G W. 

(5) A subspace V of X is called null subspace if B (v, u) = 0, for all u, v G V. 

(6) Let W be a subspace of X . The orthogonal complement ofW is the sub- 
space W 1 = {u G X | B(v, u) = 0, for all v G W}. 

(7) Let W be a subspace of X . We say that W is a non- degenerate subspace, 
relative to B, if the restriction of B to W x W is non degenerate. 

Note that the vector subspace generated by an isotropic vector u is a null 
subspace of X . 
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The next proposition states when a subspace and its orthogonal complement 
decompose the space X direct sum. 

Proposition 2.3 Let (X, B) be a generalized scalar product space of dimen- 
sion n and let W be a subspace of X . Then X = W © W ± if and only if W is 
non-degenerate. That is, the space (W, B |w) is non- degenerate, where B |w 
is the restriction of B to the subspace W. In particular, if a G X we have 
X = Ma© (Ma)" 1 if and only if B(a, a) ^ 0, where Ma denotes the subspace 
generated by a. 

For a proof see [21 Prop. 149.1]. 

In the remaider of this article (X, B) will denote a (generalized) scalar product 
space of dimension n over the field of real numbers. 

Let e = {ex, e 2 , . . . , e n } be an ordered basis of the vector space X. For each 
pair of indices i, j in {1,2,..., n} let ay := B (e iy ej) . The matrix A = [ay] 
is called the matrix of B relative to the basis e. It describes the bilinear form 
B in the following way. Let v,w in X and let x = [v] e , y = [w] e be the 
coordinate (column) vectors of v and w with respect to the basis e. Then 
B (v, w) = x t Ay. Since B is non-degenerate and symmetric, we see that A is a 
non-singular symmetric matrix. 

Proposition 2.4 Let (X, B) be a generalized scalar product space. Then there 
exist an ordered basis e* = (e{, . . . , e*) of X and nonnegative integers p and 
q, with p + q = n, such that 

(1) Bfaef) = 1 for j = 1,2,..., p. 

(2) B(e*,e*) = -1 forj=p + l,p + 2,...,p + q. 

(3) B(e*,e*)=0fort^j. 

This means that the basis e* diagonalizes the matrix associated with the 
bilinear form B. 

(4) The number of elements that satisfy B (e*, e*) = 1 is independent of the 
basis that diagonalizes the bilinear form B. 

For a proof see [2] Prop. 159.1]. 

The basis e* is called an orthonormal basis with respect to B. If q = then 
the space (X, B) is called positive definite. 

Let p and q be nonnegative integers such that n = p + q. The bilinear form 
B* on the space M n defined by 

P P+Q 
i=l i=p+l 
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where x = (xi,X2, ■ ■ ■ , x p+q ) and y = (yi, y2, ■ ■ ■ , y P + q ), is symmetric and non- 
degenerate. The generalized scalar product space (W 1 , £>*) is denoted by M. p,q . 

Any generalized scalar product space (X,B) that satisfies conditions 1,2 and 
3 of the previous proposition is isomorphic to the space M. p,q . Any such space 
(X, B) is said to have characteristic or signature (p, q), and it is usually iden- 
tified with W' q . See [2, Thm. 177.1]. 

2.1 The orthogonal group 

Among the linear operators on (X, B) the most interesting are clearly those 
that preserve the bilinear form. 

Definition 2.5 Let T : X — >■ X be a linear operator. Then T is an orthogonal 
transformation if and only if 

B (Tt>, Tw) = B (v,w) , v , w G X. 

The set of all the orthogonal transformations is a group, called the orthogonal 
group of (X, B), and denoted by O (X). 

The group O (M. p ' q ) can be considered as the set of invertible n x n matrices 
Q that satisfy Q l AQ = A, where A is the matrix associated with the bilinear 
form B with respect to the canonical basis of W° ,q . We denote O (M p '' 3 ) by 
0(p,q). 

Definition 2.6 SO (p, q) := {Q E O (p, q) \ det (Q) = 1} is the group of spe- 
cial orthogonal transformations or rotations ofM. p,q . 

If a is an invertible vector ( B (a, a) ^ 0) then the subspace (Ma)" 1 has dimen- 
sion n — 1 and it is called the hyperplane associated with a. In this case, every 
v in X has a unique representation of the form v = \a + b, with b G (Ma) and 
AgK. The next proposition shows that the linear transformation (p a : X — > X, 
defined by ip a (v) = —Xa + b is orthogonal. It is called the reflection with re- 
spect to the hyperplane (Ma)" 1 . For the sake of convenience we denote (Ma)" 1 
by H a . 

Proposition 2.7 LetW be a non- degenerate subspace of X . Define the linear 
map T : X = W © W — > X by T (v) = x — y , where v = x + y,, with xinW 
and y in W 1 - . Then T is an orthogonal transformation. 

Lemma 2.8 Let a,b be invertible vectors such that B(a,a) = B(b,b). Then 
there exists a linear map if such that (p(a) = b and (p is either the reflection 
with respect to a hyperplane or the composition of two reflections with respect 
to hyperplanes. 
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Proof: If B (a, a) — B (b, b) and a, b are invertible then a + b and a — b are 
orthogonal. This follows from B (a + b,a — b) = B (a, a) — B (6, b) . 



We have £ (a + 6, a + b) + 5 (a - b, a - b) = 4B (a, a) ^ 0. We deduce that 
either a + b is invertible, or a — b is invertible, because at least one of the 
summands in the above equation must be non-zero. 

If a — b is invertible then f a -b X ^ X is & reflection and 

Va-b (a) = ¥a-b (a - b) + ^ (a + = -1 (a - &) + ^ (a + &) = &. 

If a + 6 is invertible, then ip a +b : X —> X is a reflection and 



/ a — b a + & N 
y? 6 y? a+6 (a) = ip b ip a+b I — 1 

' a — b a + b\ 



2 

y?6 (-6) = b. 



The Cartan-Dieudonne theorem states that every orthogonal transformation 
T on an n-dimensional generalized scalar product space is the composition of 
at most n reflections with respect to hyperplanes. In the following section we 
will present our proof. 

The main difficulty to obtain the proof of the Cartan-Dieudonne theorem 
appears in the case when T(x) — x is a nonzero isotropic vector for every 
nonisotropic vector x. This case leads us to consider Artinian spaces. We 
present next some basic properties of Artinian spaces that we will use in the 
next section. 



2.2 Artinian spaces 



The simplest example of an Artinian space is the Lorentz plane R 1 ' 1 . In this 
plane the subspace generated by u = (1, 1) is a null space of dimension 1. 

An Artinian space is a generalized scalar product space of the form W' p for 
some positive integer p. Every Artinian space X has the following properties. 

(1) dim (X) = 2p is even. 

(2) X contains a null subspace of dimension p. 
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(3) If U is a maximal null subspace of X then 

(a) dim (U) = p. 

(b) If T is an element of O (X) such that T (U) = U then T is a rotation, 
that is, det(T) = 1. 

For our purposes, the main property of Artinian spaces is the following lemma 
[23 Prop. 247.1, Lemma 249.2]. 

Lemma 2.9 Let T be an element of 0(p,q) such that T(x) — x is a nonzero 
isotropic vector for every nonisotropic vector x. Then 

(1) p = q and 2p is a multiple of 4- 

(2) T is a rotation with fixed space U , where U is a maximal null subspace 
and ImiT - I) = U = Ker (T - I) . 

This result, combined with induction over the dimension of the space, is used 
in the proof of the Cartan-Dieudonne theorem presented in [21 Thm. 254.1]. In 
the next section we present an alternative proof based on an algorithm to de- 
compose a given orthogonal transformation as product of reflections. In section 
4.2 we present an explict formula for calculating the matrix representations of 
the reflections. 



3 An alternative proof of the Cartan-Dieudonne theorem 

In order to simplify the notation, from here on we write uv = B (u, v), when- 
ever B is the symmetric non-degenerate bilinear form on the space ~R p,q , and 
u and v are in W' q . The subspace generated by the vectors Ui,u 2 , ■ ■ ■ is 
denoted by (ui, u 2 , ■ ■ ■ , u^). Recall that, if a is a vector such that a 2 ^ then 
(p a denotes the reflection with respect to the hyperplane H a of all vectors that 
are orthogonal to a. 

The next lemma may be considered as a weak version of the Cartan-Dieudonne 
theorem. The proof is very similar to Uhlig's proof of the analogous result [H 
Thm. 3]. 

Lemma 3.1 Every orthogonal transformation on the space W ,q can be ex- 
pressed as the composition of at most 2n reflections with respect to hyperplanes, 
where n = p + q. 

Proof: Let T e O (p, q). Let {wi, u>2, ■ ■ ■ ,w n } be an orthogonal basis for W' q 
such that wf ^ for i = 1,2, ... ,n. Define Vj = (wj, Wj+i, . . . , w n ) for j = 
1,2,. ..,n. 

Consider the vectors T (wi) and ui\. Define the linear function ipi : ~R p ' q — > ~R p ' q 
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by 

7, if T (wx) — w 1 , 

Pi = \ <Pa, ^ T (wi) ^ w 1 and (T (u>i) - Wi) 2 ^ 0, 

^PwiPd!, if (wi) 7^ wi and (T (u>i) - u>i) 2 = 0, 
where Ci = T (wx) — w\ and di = T (w^ + w 1 . By Lemma 2.8 it is easy to see 
that ip\T (wi) = wi, and ipiT (V 2 ) C V2. 

Consider now the vectors y?iT (w 2 ) and w 2 - Define the linear function ip 2 '■ 

W' q W> q by 



V?2 



I, if (p 1 T(w 2 ) = w 2 , 

V9 C2 , if y?iT (w 2 ) ^ w 2 and (</?iT (tu 2 ) - u> 2 ) 2 7^ 0, 

^<f W2 <fd 2 , if (^2) 7^ ^2 and (</?iT (tu 2 ) - u? 2 ) 2 = 0, 



where c 2 = ip{T (u» 2 )— u» 2 and <i 2 = </?iT (w 2 )+w 2 . We know that y? 2 y?iT (u> 2 ) = 
W2, and since c 2 ,<i 2 G V 2 , we see that y? 2 (/3iT (wi) = w i- Therefore we have 
y? 2 y?iT (wi) = Wi for i = 1, 2, and y? 2 y?iT (y 3 ) C V 3 . 

Consider now (p 2 ipiT (w 3 ) and u> 3 . Define the linear function (f 3 : MP ,q — >■ M p ' 9 

by 

7, if y? 2 y?iT (w 3 ) = w 3 , 

^3 = < y? C3 , if W\T (w 3 ) ^ w 3 and (y? 2 y?iT (w 3 ) - w 3 ) 2 7^ 0, 

tfwaVda, if V?2V?iT (^3) 7^ ^3 and (y? 2 y?iT (w 3 ) - w 3 ) 2 = 0, 

where c 3 = y? 2 y?iT (w 3 ) — w 3 and d 3 = y? 2 y?iT (w 3 ) + w 3 . We know that 
y?3y?2y?i^ (w 3 ) = w 3 , and since c 3 ,d 3 G V^, we can show that ip 3 ip2fiT (w^ = 
Wi for i — 1,2, 3. Therefore we have <^ 3 <^ 2 <^iT (w^ = Wi for i = 1,2, 3, and 
(p 3 (p 2 (piT (V 4 ) C V4. 

We introduce the notation $ fe = (p k (fk~i • • -<fii lor k > I. Then, following the 
procedure used above, we can get orthogonal transformations </?i, ip 2 , ■ ■ ■ , (p n 
such that $ fc T (w^ — w { for % — 1, 2, . . . , k, and $ fc T (V k+1 ) C V fc+ i. 

The orthogonal transformations ifj are defined by 



I, if (w fc+ i) = w k+1 . 

Vc k+1 , if (tffc+i) 7^ w fc+i and ($ fe T - w k+1 f 7^ 0, 

l^wfc+iVdfc+n if ®kT (wk+i) ^ w k+1 and ($ fc T (w fe+ i) - w fc+ i) 2 = 0, 

where c k+1 = $ fc T (w k+ i) - w k+1 , and d k+1 = $ fe T (w fc +i) + w k+1 . 

Therefore we have $„T (wi) = Wi lor i = 1,2, ... ,n, and thus $ n T = I, and 
then 
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Since each ip^ is either the identity, a reflection with respect to a hyperplane, 
or the composition of two reflections with respect to hyperplanes, we see that 
T is the composition of at most 2n reflections with respect to hyperplanes and 
this completes the proof. ■ 

Reviewing the ideas used in the proof of the previous lemma, we see that one 
way to reduce the number of reflections needed to factor T is the following. 
Suppose we have found ipi,<f 2 , ■ ■ ■ ,fk, reflections with respect to hyperplanes, 
such that $ fc T (wi) — w { for % — 1, 2, . . . , £, and $ fe T (V^t-i) C V e+1 . If there 
exists Wj G {wi + i,W£ + 2, ■ ■ ■ , w n } such that $feT (wj) = Wj, or $ fe T (wj) ^ Wj 
and ($feT (wj) — Wj) 2 7^ 0, then we can reorder the elements to force j = £ + 1, 
and then, using the construction of the ipi of the previous lemma we would have 
that (fk+i is either the identity or a reflection with respect to a hyperplane. 

Lemma 2.9 tells us under what conditions it can happen that no element 
of {vj£ + i, we + 2, ■ ■ ■ ,w n } satisfies the conditions described above, and in such 
situation we can not assure that (fik+i is the identity or a reflection with respect 
to a hyperplane. 

In order to get a proof of the Cartan-Dieudonne theorem, we must find an al- 
gorithm for the construction of the orthogonal transformations ipj that avoids 
in some way reaching a situation where the hypothesis of Lemma 2.9 is satis- 
fied. This can be done by introducing an additional reflection in the way we 
describe next. 

Let T G 0(p, q) and suppose that T(x) — x is a nonzero isotropic vector for 
every nonisotropic x. Then, by Lemma 2.9 we must have p = q, n = p + q is 
a multiple of 4, and T is a rotation, that is, det(T) = 1. 

Let y be an invertible element of MP' q and let ip y be the reflection with respect 
to the hyperplane H y . Define S = (p y T. Since det(S') = —1 we see that S ^ I 
and S is not a rotation. Therefore S does not satisfy the hypothesis of Lemma 
2.9. 

Let {u>i, W2, • • • , w n } be an orthogonal basis for M. p ' q . We can reorder the ele- 
ments of the basis so that either S (wi) = W\ or S (wi) — w 1 is an invertible 
vector. Then we can find ipi that is either the identity or a reflection with 
respect to a hyperplane, and satisfies ipiS (wi) = W\ and ipiS (V2) C V2. 

Since the dimension of V2 is not a multiple of 4, the orthogonal transformation 
(fiS restricted to V 2 does not satisfy the hypothesis of Lemma 2.9. Therefore 
there exists j such that 2 < j < n and either ipiS (wj) = Wj or ipiS (wj) — Wj 
is an invertible vector. Reordering the basis of V 2 if necessary, we can suppose 
that j = 2. Then we can find ip 2 that is either the identity or a reflection 
with respect to a hyperplane and satisfies ip 2 ip\S {wi) = for i — 1,2 and 
<f2flS (V3) C V 3 . 



9 



Proceeding in the same way we can get ip^ and </?4, that are either the iden- 
tity or reflections, that satisfy {p 4 (p 3 (p 2 (piS (wi) = Wi for i = 1,2,3,4, and 
WiWzW2ViS (V 5 ) C V 5 . 

Consider now the composition where $4 = <£>4<£>3<£>2V?i- (Recall that we 
defined = ipk^k-i • • • for k > 1). There are two possible cases: 

(1) w^I for 1 = 1,2, 3, 4. 

In this case we have det($ 4 S f ) = —1 and hence $ 4 S restricted to V 5 
does not satisfy the hypothesis of Lemma 2.9. Therefore we can find 
orthogonal transformations tpj, for j = 5, 6, 7, 8 such that &sS (wi) = Wi 
for i = 1, 2, . . . , 8 and &$S (Vg) C Vg. Then, for the composition <J> 8 5' we 
have again the same two possible cases that we had for but now 
considering the maps pi for i — 5, 6, 7, 8. 

(2) tfi — I for at least one i, with 1 < % < 4. 

In this case dim (V5) is a multiple of 4 and it is possible that det = 

1. 

If restricted to the space V 5 does not satisfy the hypothesis of 
Lemma 2.9 then we can find <^ 5 , which is either the identity or a reflection 
such that $ b S{wi) — Wi for % — 1, 2, . . . , 5 and & 5 S (V 6 ) C \/ 6 . Notice 
that, since at least one of the 93^ is the identity, for 1 < % < 4, the number 
of reflections in the composition $ 5 is at most equal to 4. 

If $ 4 S' restricted to the space V 5 satisfies the hypothesis of Lemma 
2.9 then we take an invertible vector z in V 5 and form the composi- 
tion (fiz^^S. This map restricted to V 5 can not satisfy the hypothesis of 
Lemma 2.9, and consequently, we can find tp^, which is either the iden- 
tity or a reflection and satisfies ip^ip z ^ 4 S (wi) — Wi for i — 1, 2, . . . , 5 
and psPz^aS (Vq) C Vq. Since at least one of the (fi is the identity, for 
1 < i < 4, the number of reflections in the composition <y9 5 <^<I>4 is at most 
equal to 5. 

Applying the procedure described above we see that for each £ such that 
1 < £ < n, we can find reflections with respect to hyperplanes ipi, ip2, ■ ■ ■ , (p s 
such that $ S S (vi) — i>j for 1 < i < £, where {vi, . . . ,v n } is a reordering 
of the orthogonal basis {w±, . . . ,w n }, and s < £. This last inequality is very 
important. 

In particular, for £ = n we get <fr s S = I, with s < n. We claim that the case 
s = n is not possible. If s — n then det (<& s ) = (— l) n = 1, because n is a 
multiple of 4. On the other hand, det (<& s ) det(S') = det(J) = 1. But we know 
that det(S') = — 1. Therefore s = n is not possible and we conclude that s < n. 

Since S = <p y T ', we have <& s (p y T = I and therefore T is the composition of at 
most n reflections with respect to hyperplanes. 

We have proved the following result. 
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Lemma 3.2 LetT be an element ofO(p,q). IfT (x)—x is a nonzero isotropic 
vector for every nonisotropic vector x then T is the composition of at most 
p + q reflections with respect to hyperplanes . 

Lemma 3.3 Let T be an element ofO(p, q) . If there exists a basis {w±, . . . , w p+q }, 
where all the elements are nonisotropic, such that 

T (wi) - Wi 

is a nonzero isotropic vector for i — 1, . . . ,p + q, then T is the composition of 
at most p + q reflections with respect to hyperplanes. 

Proof: We can proceed as follows. 

Step 1 In each step we have to deal with three possible cases. In particular, 
here we have: 

(1) There exists a nonzero and nonisotropic element v 1 G W p ' q such that 
T(v 1 )-v 1 = 0. 

(2) There exists a nonzero and nonisotropic element v\ of the basis such that 
T (vi) — v i 7^ and (T (vi) — t>i) 2 ^ (i.e. T (vi) — Vi is nonisotropic). 

(3) For every nonzero and nonisotropic x G W' q we have that T (x) — x ^ 
is isotropic. 

In case 3), from Lemma 3.2 we obtain that T is the composition of at 
most p + q reflections with respect to hyperplanes. 

In cases 1) or 2), we can find ipi, which is either the identity or a reflection, 
that satisfies 

<PiT(ui) = vi 

and also ipiT (Wi) = Wi, where W 1 = (wi)- 1 . 
Step 2 Now, for the orthogonal transformation ipiT restricted to the space 
Wi, we have 

(1) There exists a nonzero and nonisotropic element i>2 G W\ such that 
<PiT(v 2 ) - v 2 = 0. 

(2) There exists a nonzero and nonisotropic element v 2 G W\ such that 
PiT (v 2 ) — v 2 7^ is nonisotropic. 

(3) For each nonzero and nonisotropic x G M. p ' q we have that (p\T (x) —x ^ 
is isotropic. 

In case 3), we have that (p{F = S and thus, from Lemma 3.2, S is the 
composition of at most p + q — 1 reflections with respect to hyperplanes. 

In cases 1) or 2), we can find tp 2 , which is either the identity or a reflection, 
that satisfies 

(p2<PlT(Vi) = Vi, 

where % can be chosen as 1 or 2. Also, it is fulfilled that ^p 2 ^iT (W 2 ) = W 2 , 
where W 2 = (wi,W2) ± - 

By following these steps, we can end up with one of the two following situa- 
tions: 
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A We can find an orthogonal set of nonisotropic elements {vi,v 2 , • • • ,v p+q } 
and a finite sequence of linear transformations ip±, tp 2 , • • • , Pp+q, such that: 

• (fi is either a reflection or the identity, for % — 1, 2, . . . ,p + q. 

• tpitpi-x ■ ■ ■ (vi) = v u for i = 1, 2, . . . , I and I = 1, 2, . . . p + q. 

• ipiipi-i ■ ■ ■ (P2¥\T (Wi) = Wi, for I = 1,2, .. .p+q, where Wi = (vi, v 2 , ■ ■ ■ , Vi) 
for i — 1, 2, . . .p + q. 

B We can find an orthogonal set of nonisotropic elements {i>i, i>2, . . . , Vk} and 
a finite sequence of linear transformations (p±, <p 2 , ■ ■ ■ , fk, where k < p + q, 
such that: 

• ifii is either a reflection or the identity, for i — 1, 2, . . . , k. 

• tpitpi-i ■ ■ ■ y?2V?iT (v^ = Vi, for i = 1, 2, . . . , / and I = 1, 2, ... A;. 

• ^/^/-i ' ' • V^i^ (Wj) = Wj, for / = 1, 2, . . . k, where = (f 1, v 2 , . . . , Wj)" 1 
for i = 1, 2, . . . A;. 

• Consider the orthogonal transformation y^y^-i ■ • • (p 2 (piT, restricted to 
Wk- Then we have that for every nonzero and nonisotropic x G it is 
fulfilled that (pk<pi^i ■ ■ ■ (p 2 (fiT (x) — x is nonzero and isotropic. 

In case A, we have that 

Wf sVVh?-i ' ' ' ^PiT = I. 
Thus T is the composition of at most p + q reflections. 
In case B, we have that 

<£W-i ' ■■W2 ( fiT = S, 

restricted to Wk (whose dimension is p + q — k). From Lemma 3.2 we conclude 
that S is the composition of at most p + q — k reflections. 

In summary, T is the composition of at most p + q reflections. ■ 

Now we are ready to prove the Cartan-Dieudonne theorem. 

Theorem 3.4 Let T be an element of 0(p,q). Then T is the composition of 
at most p + q reflections with respect to hyperplanes. 

Proof: Let n = p + q and let {w±,w 2 , . . . , w n } be an orthogonal basis for W' q 
with wf 7^ for 1 < % < n. 

If there exists an element Wj of the basis such that T (wj) = Wj or T (wj) 7^ Wj 
and (T (wj) — Wj) 2 7^ 0, then, reordering the basis (keeping the notation Wi 
for the basic elements), we can get j = 1. Thus we can find ip±, which is either 
the identity or a reflection, that satisfies (p\T (wi) = wi, and Lp{F (V 2 ) C V 2 . 

If there exists an element Wj G {^2, ^3, • • • , w n } such that ipiT (wj) = Wj or 
ipiT (wj) 7^ Wj and (tpiT (wj) — Wj) 2 7^ 0, reordering the basis {w 2 , w 3 , . . . , w n } 
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we can assume that j = 2, and then we can find </?2> that is either the identity 
or a reflection with respect to a hyperplane, that satisfies ip2<piT (u>i) = Wi for 
i = l,2, and (p 2 ipiT (V 3 ) C V 3 . 

This process can be repeated until we reach either of two possible cases. In 
the first case we can find tpi, for 1 < i < n, such that $ n T = /, and therefore 
T is the composition of at most n reflections with respect to hyperplanes. 

In the second case we find (fi, for 1 < i < k, where 1 < k < n, that satisfy 
§ k T(wi) = Wi for 1 < i < k, <f> k T (Vk+i) Q V k+1 , and ® k T(wj) - w j is a 
nonzero isotropic vector for every para k < j < n. Applying Lemma 3.3 to 
the orthogonal transformation $fcT, on the space Vk+i, we obtain reflections 
with respect to hyperplanes Tj, for 1 < i < s, with s < n — k, such that 
t~ s t~s-i ■ • -Ti$feT is the identity on W ,q . Since s + k < n, we conclude that, 
also in this case, T is the composition of at most n reflections with respect to 
hyperplanes. ■ 



4 Clifford algebras 

In this section we present some definitions and basic results about Clifford 
algebras associated with a generalized scalar product space (X, B) of signature 
(p, q). We use the notation of [7J Ch. I]. We also describe how the Clifford 
algebra structure can be used to deal with orthogonal transformations, and in 
particular, with reflections with respect to hyperplanes. 

Definition 4.1 Let (X,B) be a generalized scalar product space of dimension 
n and let A be a real associative algebra with identity 1 such that 

(CI) A contains copies o/R and of X as linear subspaces. 

(C2) For all v G X we have v 2 = B (v, v) . 

(C3) A is generated as a ring by the copies ofM. and of X , or equivalently, as 
a real algebra by {1} and X . 

Then A is called a real Clifford algebra for (X, B) and it is denoted by A = 
C(X). 

Note that (C2) links the multiplication in the algebra with the bilinear form 
on X. 
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4-1 Bases for Clifford algebras 

Here we describe a basis for C (X) in terms of an orthonormal basis for X. 
Let (X, B) be a generalized scalar product space of signature (p, q) and let 
e = {e±, e 2 , . . . , e„} be an orthonormal basis for X. Then, by (C2) we have 

1,2,. ..,p 

p + l,p + 2,...,p + q, 
and 

- (v w + wv) = B (v,w) , u,v G A\ 
Define iV = {l,2,...,n}. Let ft, ft, • • • , ft be distinct elements of N. Then 

e /3i e /3 2 • • - e /3 3 = (-l) a e ai e a2 • • -e aa , 

where (a±, a 2 , • • • , « s ) is the permutation of 02, ■ ■ ■ , f3 s ) that satisfies «i < 
0:2 < • • • < «n, and o" is the number of transpositions of the permutation that 
sends (a ± , a 2 , ■ ■ ■ , a s ) to (ft, ft, • • • , ft)- 

Since ef = ±1, it is easy to see that, up to a change of sign, every product of 
basic elements, possibly with repeated factors, can be reduced to a product 
of at most n factors with indices in increasing order. This implies that every 
element of C (X) can be written in the form 

A 

where the sum runs over all the subsets A = {a±, a 2 , ■ ■ ■ , a n } of iV (with 
«i < a 2 < ■ ■ ■ < a„) and is a real coefficient. In order to simplify the 
notation we define = e ai e Q2 • • • e as , and we put eg = 1. Therefore the 
collection of elements of the form e^, where A is a subset of N, is a generating 
set for C (X) and consequently, dim [C (X)) < 2 n . 

Recall that a vector s e X is called invertible if B (s, s) 7^ 0. This condition 
is equivalent to s 2 = B(s,s) 7^ 0. Therefore, the element (1/B (s, s)) s is the 
inverse of s with respect to the multiplication in the algebra C(X), and we 
can write s" 1 = s/s 2 = s/B(s,s). 

Theorem 4.2 Let C (X) be a Clifford algebra for an n- dimensional general- 
ized scalar product space (X, B) with signature (p, q) and let 

e = {e±, . . . , e p , e p+ i, . . . , e p+q } be an orthonormal basis for (X, B). Then 



e 2 = B(e h d) = 
and it is easy to show that 
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(i) If n is even then dimC (X) = 2 n and {ca '■ A C N} is a basis for C {X). 

(ii) If n is odd and ^ M, then dimC (X) = 2 n and {ca '■ A C A^} is a basis 
forC(X). 

(Hi) If n is odd and G R , then = ±1 and p — q = 1 mod (4). In this 

case dimC (X) = 2 n ~ 1 and {e^ : A C. N, #A even } is a basis for C (X). 

Theorem 4.3 Let C (X) be a Clifford algebra for the generalized scalar prod- 
uct space (X,B) and let W be a non-degenerate subspace of X. Then the 
subalgebra of C (X) generated by W is a Clifford algebra for W . 

The proof of these results can be found in [7J, [8], and [9]. 

A result of Chevalley states that every n-dimensional generalized scalar prod- 
uct space (X, B) has a Clifford algebra of dimension 2 n (see [7]). If (X, B) has 
signature (p, q) then it is isomorphic to W ,g . The Clifford algebra of dimension 
2 n of W'i is denoted by R Pi5 . 

The algebra R M is a graded vector space. Let s be an integer such that 
< s < n. Let denote the subspace generated by {ca '■ A C N, j^A = s}. 

It has dimension ( "J . The elements of M.^ q are called s- vectors. Note that 

n 
s=l 

The space of 0- vectors is generated by {1}, and its elements are also called 
scalars. The space M p can be identified with M p ' g . The space is generated 
by {ejv}. The element ex is called pseudoscalar. 

Every element a G IR p>g can be expressed in the form 

a = ^A&Ai Xa e 
ACiV 

and also as 

n 

a = J2a r , a r G K£ j9 . 

r=0 

Definition 4.4 Let a = I]" =0 a r; where a r G R» ? - If a r = for r > t and 

a t 0, then we say that the grade of the multivector a is t, and write gr[a) = t. 

Notice that if A G M then gr(Aa) = gr(a). If si, . . . , are non-zero vectors in 
M. p,q , then gr (si • • • Sk) = I < k and dim(si, . . . , S&) = I. 

The multiplication in the algebra IR Pjg can be used to represent reflections with 
respect to a hyperplane in the space ~R p ' q as follows. Let s be an invertible 
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vector in W> q . Define the map (f s : W' q ->■ W> q by 

^ s ( x ) = -sxs' 1 , x e R p ' q . 



In order to show that tp s (x) is a vector in MP ,q for all a; G ~R p ' q we first compute 
sxs and get 



sxs = — (sx + xs + (sx — xs)) s 
= - (2B (x, s) + sx — xs) s 
= B (x,s) s + -sxs — -xs 2 . 

The last equation yields sxs = 2B (x, s) s — xs 2 and hence 

-i 2£(x,s) 
— sxs = S + X, (1) 

s 2 

which is clearly a vector in M. p,q . 

Lemma 4.5 Let s be an invertible vector in W ,q . Then the linear map ip s : 
MP' q — > MP' q is an orthogonal transformation. Furthermore, it is the reflection 
with respect to the hyperplane H s = {x e M, p,q : £> (x, s) = 0}. 

Proof: For x, y G W' q we have 



i3(y? s (x) , Lp s (y)) =B (-sxs 1 ,-sys ^ 

= \ ((- sxs ^ 1 ) (-sy^ 1 ) + (sys' 1 ) (-sxs -1 )) 

= 7^ (xy + yx) s -1 
= B(x,y) . 

If x e H s we have (x) = — sxs -1 = — (2B (x, s) / s 2 ) s + x = x. 

If x = As where A is a real number then ip s (x) = — s (As) s _1 = —As = — x. 
Therefore cp s is the reflection with respect to H s . ■ 

It is easy to verify that ip s satisfies 

(1) ip s = ipx s for every nonzero real A. 

(2) y?; 1 = yv 

(3) If si, s 2 , ■ ■ ■ , Sfc are invertible elements of W ,q then 



<Psi<Pn ■■■ ( fs k (x) = S!S 2 • • • Sfc.xSfcXSfc 1 S k ^ 1 ■■■S 2 1 S 1 1 
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Using the Clifford algebra multiplication the Cartan-Dieudonne theorem reads 
as follows. 

Theorem 4.6 Let T be an orthogonal transformation on M. p,q . Then there 
exist invertible elements Sx, s 2 , ■ ■ ■ , in ~R p ' q , with k < p + q, such that 

T (x) = (-l) fe sxs 2 ■ ■ ■ Sk-xSkXs^s^x ■ • ■ s^s^ 1 , x e W ,q . 



If T is an orthogonal transformation on IR P,<7 and there exist invertible multi- 
vectors A, B in the Clifford algebra such that 



T (x) = ±AxA~ 1 = ±BxB-\ x 6 R M , 

then A = XB for some real A. The proof of this result can be found in [9]. We 
will illustrate this fact in the examples of the next section. 



4-2 A matrix representation of ip s 



Having in mind that, in many applications, matrix representations of reflec- 
tions are useful, we introduce another algebraic expression for ip s . 

Define A = then 



s = Xx + y, 

where y is orthogonal to x, that is 

xy = —yx. 

Thus 



— sxs = — - (Ax + y) x (Xx + y) , 
s 2 

= ~ ((AV - y 2 ) x + 2Ai- 2 y) , 

=-?(( l ¥ : -» , ) ,+i »^)')' (2) 

This last result turns out to be useful to find a matrix representaction of 
ip s with respect to a given orthogonal basis B = {ui\, w 2 , ■ ■ ■ , w p+q }, where 
wf 7^ for i = 1, 2, . . . ,p + q. Indeed, consider s = Y%=x a i w i, then 
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p+q 



s = a k w k + a > n ' 
1=1 

i^k 



t=l 



Define 



2 W k + V 5 «J». 



P+9 P+<3 



Vk = 2^ OLiWi = > g Wi 

i=l i=l 



Using ()2]) to get 



(( 



f (B(s,w k )f 
{ < 

(B(s,w k )) 2 



yl)w k + (2B (s,w k ))y k 



( 



Wi 



p jZ{B{s, Wl )) 
. »=i 



»'7 



w k + 2B(s,w k )y k 



J J 



J 



From this result we obtain that the k-th column of the matrix A s = [(p s ] B is 
given by 



[ifs (W k )] 



B 



Equivalent ly 



2B(s,w k ) 



2B(s,w k )^ 



(B(s,w k )) 2 _ P X q (B(s,w t )f 
"A i= 1 



2B(s,w k ) 



B{s,w p+q ) 

2 

W 

v+q 



(fc-th row) 



(^)« 



(S(s,^)) 2 



for I = j. 



(3) 
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5 Examples 



Let T be the orthogonal transformation on the space R 2 ' 3 represented by the 
matrix 



1 5 4 3 
-5 1 3-4 
4 3 15 
3-4-5 1 
-1 



with respect to the canonical basis E = {ei, e 2 , e 3 , e±, e$} of IR 2,3 . Then we have 
T (ei) - ei ^ 0, (T (e^) - e^) 2 = for i = 1, 2, 3, 4, and T (e 5 ) - e 5 = -2e 5 . 
We can take c\ = e 5 . 

It is easy to see that (p Cl T restricted to (ei, e 2 , e 3 , e±) satisfies the condi- 
tions of Lemma 2.9. By Lemma 2.8 and the proof of Lemma 3.1 we obtain 
c 2 = faT (ei) + e\ = 2ei — 5e 2 + 4e 3 + 3e 4 and c 3 = e\. Then we have 
{p C3 ip C2 {p Cl T (e^ = ei, for i — 1,5. Notice that ip C2 ip Cl T does not satisfy the 
hypothesis of Lemma 2.9 on (ei, e 2 , e 3 , e^). 

Now we have (ip C3 ip C2 ip cl T (e») — ei) 2 ^ for i = 2, 3, 4. We can take 



25 23 

C4 = <p<sip C2 <p Cl T (e 2 ) - e 2 = ye 2 - 7e 3 - ye 4 . 



Then we have ip Ci ip Ci ip C2 ip Cl T (e^) = for i = 1, 2, 5, and 



((p Ci (p Cz ip C2 (p Cl T (e^ - eif ^ for i = 3, 4. 

We can take c 5 = (p C4 tp c3 (p C2 (p Cl T (e 3 ) - e 3 = ^e 3 + ff e 4 - 



It is easy to verify that ip C5 ip C4 ip C3 ip C2 ip Cl T (ei) = ei for i = 1,2,3,4,5, that is 
T = (p cl <p ca <Pca<Pc i <Pc S - 



Let us denote by Aj the matrix representation with respect to the canonical 
basis E of <p c ., for 1 < j < 5. Using formula we obtain: 
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1 













1 
















-1 


1 
















27 
2 


7 


23 
2 







10 





7 

25 


24 
25 





, A 4 = 





-7 


73 
25 


161 
25 





, ^3 = 


10 





24 
25 


7 
25 










23 
2 


161 

25 


479 
50 







1 











1 
















1 




1 



-1 


5 


4 


3 







10 





5 


23 
2 


-10 


15 

2 







10 





-4 


10 


9 


6 





, A 1 = 


10 





-3 


15 
2 


6 


n 

2 







1 

















1 







-1 



A direct computation yields T E = A 1 A 2 A 3 A 4 A 5 . 

Computing the Clifford product C1C2C3C4C5 we obtain a linear combination of 
s-vectors where s < 3. This suggests that C1C2C3C4C5 may be equal to the 
Clifford product of 3 vectors and that T could be factored as the product 
of 3 reflections. We present next another way to factor T that confirms this 
conjecture. 

We apply the factorization algorithm to T, but now using the orthogonal basis 
W = {wi,W2,ws,W4,W5}, where 

wi = e 3 + e 4 -e 5 , w 2 = e 1 + e 2 ,w 3 = e 1 + e 4 + 2e 5 ,w 4 = e 3 -e 4 and w 5 = ei-e 2 . 
We obtain (T (wi) — Wi) 2 7^ 0, for 1 < % < 5. Therefore we can take 



d l =T{w 1 ) - W!, 

ip dl T (wi) -Wi^O and ((p dl T (wi) - Wif ^ 0, for i = 2, 3, 4, 5, 
d2 = <Pd 1 T (w 2 ) -w 2 , 
(p d2 (p dl T (wi) -Wi^O and (ip d2 ip dl T (1^) - wif ^ 0, for i = 3, 4, 5, 
d3 = ¥d 2 ¥d 1 T {w 3 ) -w 3 , 
= ip d3 (p d2 ip dl T (wi) - w u for i = 4, 5. 

and therefore T = <f dl <f d2 ip d3 . 

Let .Bj denote the matrix representation of ip di with respect to the canonical 
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basis E, for i — 1,2, 3. We have 



51 


7 


35 


35 


2 


2 


2 


2 


7 


3 


5 


5 


2 


2 


2 


2 


35 


5 


23 


25 


2 


2 


2 


2 


35 


5 


25 


23 


2 


2 


2 


2 


7 


-1 


-5 


5 



di = 7e 1 - e 2 + 5e 3 - 5e 4 + 2e 5 , 



347 


104 


286 


208 


26 


9 




9 


9 


9 


3 


104 


41 


88 


64 


8 




9 


9 


9 


9 


3 


286 


88 


233 


176 


22 


9 




9 


9 


9 


3 


208 


64 


176 


119 


16 




9 


9 


9 


9 


3 


26 


8 


22 


16 


-1 


3 




3 


3 


3 






43 


25 


35 


5 


5 




18 


18 


18 


18 


3 




25 


43 


35 


5 


5 




18 


18 


18 


18 


3 




35 


35 


31 


7 


7 




18 


18 


18 


18 


3 




5 


5 


7 


17 


1 




18 


18 


18 


18 


3 




5 


5 


7 


1 


-1 




3 


3 


3 


3 





d 2 = 26ei — 8e 2 + 22e 3 — 16e 4 + 6e 5 , 



d 3 = -5ei + 5e 2 - 7e 3 + e 4 - 6e E 



It is easy to verify that T# = BiB 2 B 3 , and notice that d\d 2 d 3 = §c\c 2 c 3 Ct±c^. 
The matrix representation of T with respect to the basis W is 



w 



I 9 i _10 8 

3 3 3 3 



3 13 4-5 

2 -i _1 _5 4 

3 3 3 3 

5 4 5 1 -3 

4 5 4 -3 1 



Note that this second factorization avoids the Artinian case, that is, the situ- 
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ation where the hypothesis of Lemma 2.9 holds. 



As another application of formula (|3]), we can find the matrix representation 
of each ip di , for i = 1,2,3, but now with respect to the orthogonal basis W. 
Using the above values of di, d 2 and d 3 , respectively, we obtain: 



1 -A 

- —2 


A 
1 


i n 

JLU 


O 
o 


3 Z 


3 


T 


3 


3 1U 


— D 


— 15 


1 o 


i 2 


1 

± 


1 n 

1U 


Q 
o 




~~ 3 


3 


3 


5 15 


-10 


-24 


20 


4 12 


-8 


-20 


17 


1 











10 


-9 


-19 


17 


3 


-2 


19 
3 


17 
3 


19 


-19 


352 
9 


323 
9 


17 


-17 


323 
9 


298 
9 



and 





1 













1 











c 3 = 





-2 


4 
3 


5 
3 








-4 

-5 


7 
9 

20 
9 


20 
9 

34 
9 



A direct computation yields T w = CiC 2 C 3 . It should be noticed that the 
orthogonal basis W is not ordered according to the signature of IR 2 ' 3 . 



6 Final remarks and conclusions 

The factorization of an orthogonal transformation is not unique, as we have 
shown in our example. Actually, the number of reflections that factorizes an 
orthogonal transformation is not unique either. An interesting question related 
to this last point is to determine the minimum number of reflections by hy- 
perplanes required to factorize a given orthogonal transformation. An answer 
to this question using matrices has been given in [10], and it was translated 
to the language of orthogonal transformations in (21 Pgs 260,261]. 
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In our examples (Section [5]) we use the grade of a multivector to find the 
minimum number of reflections by hyperplanes required to factorize Tg. We 
can formalize this procedure as follows: 

Lemma 6.1 Let T be an orthogonal transformation on M. p,g . Assume that 
there exist invertible elements s\, . . . , Sk in MP' q , with k < p + q, such that 

T (x) = (-l) fc sis 2 ■ ■ ■ Sk-xSkXs^s^ ■ ■ ■ s 2 1 s^ 1 , x G W> q . 

If gr (si ■ ■ ■ Sk-iSk) = t < k, then T cannot be factored into less than t reflec- 
tions with respect to hyperplanes. 

Proof: Suppose that there exist s[, . . . , s[, with / < t, such that 

T (x) = (-l) fc 4 • • ■ sl^xis'^isl,)- 1 ■ ■ ■ (s[)-\ x E BM. 

We have that gr (^s[ ■ ■ ■ s' l _ 1 s' l ^j < I < t. But (^s[ ■ ■ ■ s'^s'^j = A (si • • • Sk-iSk), 

where A G R and therefore gr (s^ ■ ■ ■ sj^sj) = gr (s% ■ ■ ■ Sk-iSk) = t, which is a 
contradiction. ■ 

With this result, we can state the following theorem 

Theorem 6.2 Let T be an orthogonal transformation on the space R p ' 9 . // 
there exist invertible elements Si, s 2 , ■ ■ ■ , Sk G R p,9 ; such that: 

• T(x) = (-l) fc s x ■ ■ ■ Sk-iSkXs^s^ ■ ■ ■ s 2 1 s± 1 , x G W' q , 

• gr(s 1 ,...,s k „ 1 ,s k ) =t, and 

• KeriT — I) is non-degenerate, 

then T can be factored into t reflections with respect to hyperplanes and, more- 
over, dim(Ker(T — I)) 1 = t. 

Proof: By hypothesis T = <p sl ---<p Sh and gr (si, . . . , s fc _i, s k ) = t. By con- 
sidering the space V% = (s±, ■ ■ ■ ,Sk), where dimVi = t, it is easy to show 
that C Ker (T — I). If Ker (T — /) is non-degenerate. Then we can find 
an orthogonal basis B = {wi, . . . ,Wj, Wj + i, . . . , Wj + {\, where j + I = n and 
Ker (T — I) = (wi, . . . , Wj-i, Wj) and applying the algorithm proposed in this 
paper we get tpi = I, for i = 1,2, ... ,j. From the application of the Cartan- 
Dieudonne theorem to T\v J+1 , where Vj+i = (wj+i, . . . ,Wj + i), we obtain that 
T is the composition of at most I reflections by hyperplanes. From Lemma 6.1 
we know that t < I. 

Now, suppose that there exist u\,...,u m G (Ker (T — I)) ± , such that T = 
Pui^'Vum- If g r (^i""" M m) < I then by considering T\ Vj+1 we can find a 
nonzero vector u G Vj + i orthogonal to the set {u\, . . . u m }. But in such case 
u G Ker (T — I) D (Ker (T — I)) , which is a contradiction since Ker (T — J) 
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is non-degenerate. Therefore T is the composition of I reflections through 
hyperplanes, that is, there exist U\,...,ui such that T = <£> U1 • • • and 
gr(ui---ui) = I and, moreover, t = gr (sis 2 • • • Sfc) = gr (uiu 2 • • • ui) = I = 
dim(Ker(T-7)) ± ■ 

It should be remarked that the previous Theorem relates the grade of a mul- 
tivector with the Cartan-Diedudonne-Scherk theorem. It remains however, to 
propose an algorithm capable of finding explicitly the minimum number of 
reflections required to decompose a given orthogonal transformation. 
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